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We present a generalization of the electron localization function (ELF) that can be used to analyze 
time-dependent processes. The time-dependent ELF allows the time-resolved observation of the 
formation, the modulation, and the breaking of chemical bonds, and can thus provide a visual 
understanding of complex reactions involving the dynamics of excited electrons. We illustrate the 
usefulness of the time-dependent ELF by two examples: the tv-tv* transition induced by a laser field, 
and the destruction of bonds and formation of lone pairs in a scattering process. 
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The advent of ultrashort laser sources with pulse dura- 
tions on the order of 10-100 femtoseconds [l| has paved 
the way to the analysis and control of chemical reac- 
tions 2]: By means of pump-probe spectroscopies with 
femtosecond laser pulses |2| one can follow, in the time 
domain, the nuclear motion which typically evolves on 
the picosecond time scale. One of the most important 
recent achievements has been the experimental realiza- 
tion of attosecond pulses 0. These are produced by co- 
herently superimposing high harmonics generated by a 
strong infrared laser pulse interacting with atoms. With 
this light source available, pump-probe spectroscopies us- 
ing attosecond pulses allow the temporal resolution of 
the electronic motion which is governed by the femtosec- 
ond time scale. Questions like "How does an electron 
travel from the highest occupied to the lowest unoccupied 
molecular orbital when excited by a laser?" may soon 
become experimentally accessible. This establishes the 
need for theoretical tools to analyze and interpret such 
data. Theoretical tools of this type will be developed in 
this article. We are mainly concerned with bonds, and 
how they break, form, or change during a time-dependent 
process. 

The intuitive concept of a chemical bond is very simple 
and elegant: an electron pair shared between neighbor- 
ing atoms that provides the necessary attraction to bind 
the molecule. However, it turns out to be very difficult 
to define exactly what a bond is, or even to visualize 
it. The one-electron molecular orbit als that stem from 
density-functional theory or Hartree-Fock usually have 
contributions from several atoms and do not represent a 
unique bond. The electronic density, on the other hand, 
does not easily reveal important features like lone pairs. 
The electron localization function (ELF) is a function 
crafted to bring into evidence the subtle bonding proper- 
ties of an electronic system. It was originally applied to 
ground- state systems, in the study of atomic shells and 
covalent bonds 5]. Soon after, it was realized that the 
ELF could be also used to analyze lone pairs, hydrogen 
bonds |(|, ionic and metallic bonds 0, etc. The systems 



studied include atoms 0, molecules 0, surfaces and 
solids H, 0. It is also possible to establish a rigor- 
ous topological classification of chemical bonds using the 
ELF [l£|. Furthermore, the ELF has the advantage of 
being fairly insensitive to the method used to calculate 
the wave functions of the system. In fact, Hartree-Fock, 
density-functional theory, or even simple approaches such 
as extended Hiickel methods, yield qualitatively similar 
ELFs 7]. Approximate electron localization functions 
have also been obtained from experimental electron den- 
sities measured with x-ray 0. 

Up to now the ELF has only been used to study sys- 
tems in their ground state. Clearly, an extension of the 
ELF to time-dependent processes appears quite desir- 
able. Such an extension would allow the time- resolved 
observation of the formation, the modulation, and the 
breaking of chemical bonds, and thus provide a visual 
understanding of complex reactions involving the dynam- 
ics of excited electrons. That is exactly the goal of this 
communication. 

The derivation of a time-dependent ELF (TDELF) fol- 
lows closely the derivation of Becke and Edgecombe of the 
ground-state ELF 5]. Our objective is to find a func- 
tion, D a (r,t), that measures the probability of finding 
one electron in the near vicinity of a reference like-spin 
electron at position r and time t. If this probability is 
high then the reference electron must be delocalized. On 
the other hand, a covalent bond consists of a localized 
electron pair (of opposite spins) between two neighbor- 
ing atoms. The probability of finding a like-spin electron 
close to a reference electron in this region must then be 
very low. For a determinantal many-body wave function, 
the function D a (r,i) is given by 



4 n a (r,t) n a (r,t) 



(1) 



where a denotes the spin, n a the spin density, j a the 
absolute value of the current density, and 



r a (r,t) = ^|V(/^(r,£)| 2 . 



(2) 



i=l 
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This last expression represents the kinetic-energy density 
of a system of N a electrons, described by the one-particle 
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FIG. 1: (Color) Snapshots of the time-dependent ELF for the excitation of acetylene by a 17.15 eV laser pulse [l^. The pulse 
had a total length of 7 fs, an intensity of 1.2 x 10 14 Wcm -2 , and was polarized along the molecular axis. Ionization and the 
transition from the bonding tt state to the antibonding 7r* state are clearly visible. 



orbitals (fi a . These orbitals can be obtained, e.g., from 
time-dependent density-functional theory or from a time- 
dependent Hartree-Fock calculation. 

Equation (QJ is similar to the expression for the 
ground-state D a (r) p. The main difference is the ad- 
ditional term proportional to j% This term natu- 
rally arises when the analysis of Ref. 5] is carried out 
without assuming real- valued orbitals, i.e., without as- 
signing vanishing orbital currents. This was pointed out 
by Dobson [ll[ in his evaluation of the Fermi-hole curva- 
ture of a static current-carrying single-determinant state. 
It is easy to see that the derivation of Ref. [ll| applies 
equally well to time-dependent determinantal wave func- 
tions, thus leading to Eq. (^J. 

The function D a (r,t) is always > 0, but it is not 
bounded from above. As usual, we define as an alter- 
native measure of localization 



l+[D a (r,t)/D0(r,t)} 2: 
with the definition D%(r,t) = r™ G (n a (r, t)) , where 



(3) 



HEG 



K) = -(6tt 2 ) 



2/3 n 5/3 



(4) 



is the kinetic-energy density of a homogeneous electron 
gas of density n a . Using this definition, /elf is dimen- 
sionless and lies between zero and one. A value of 1, i.e., 
D a (r,t) approaching zero, corresponds to high localiza- 
tion. 

We illustrate the usefulness of the time-dependent ELF 
by two examples: (i) the excitation of acetylene by a 
strong laser pulse (Fig. [Q, and (ii) the scattering of a 



high-energy proton from the ethene molecule (Fig. [2J. 
The figures depict a slab of the ELF passing through the 
plane of the molecules, to which we superimposed an iso- 
surface (Fig. [Q or contour lines (Fig. |2J) at /elf = 0.8. 
Movies of the time-dependent ELF and of the corre- 
sponding time-dependent density can be found in our 
website [l3j or as supplementary material to this article 
|l2|. All calculations were performed in the framework 
of time-dependent density- functional theory [3, , us- 
ing a real-space, real-time code |l6j |. and employing the 
adiabatic local-density approximation ^| for the time- 
dependent exchange-correlation potential. This approxi- 
mation is expected to overestimate ionization due to spu- 
rious self-interactions. This problem can be overcome by 
the use of explicitly orbital-dependent functionals within 
the time-dependent optimized effective potential method 

m. 

In the beginning of the simulation (i) the acetylene 
molecule is in its ground state [Fig. ^a)]. At this mo- 
ment, the ELF exhibits three major features: a torus 
between the carbon atoms — typical of triple bonds — and 
the two characteristic blobs around the hydrogens. As 
the intensity of the laser increases, the system starts to 
oscillate, and then to ionize [Figs, ^b) and^c)]. Note 
that the ionized charge leaves the system in fairly local- 
ized packets [the blob on the left in panel (b), and on 
the right in panel (c)], that then spread with time. The 
central torus then starts to widen until the moment it 
breaks into two separate tori, each around one carbon 
atom [Fig. [Je)]. We interpret this finding as a transition 
from the tt bonding state into the 7r* nonbonding state. 
The system then remains in this excited state for some 
time. We emphasize that our calculation corresponds to 
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FIG. 2: (Color) Snapshots of the time-dependent ELF for the scattering of a fast, nonrelativistic proton (E^i n ~ 2 keV) 
by ethene |l2(. The molecule breaks into several pieces. During this fragmentation process, the breaking of bonds and the 
subsequent creation of several lone pairs become clearly visible in the time-dependent ELF. The legend is the same as in Fig.[l] 



one specific orientation of the molecule, namely the ori- 
entation where the polarization vector of the laser field 
is parallel to the molecular axis. A detailed analysis of 
the dependence of photoionization on the molecular ori- 
entation will be presented elsewhere pj| . 

To assess the importance of the j 2 /n term, we com- 
pare, in Fig. 03 a snapshot of the proper time-dependent 
ELF [given by Eq. Q] with the traditional ground-state 
expression (where the j 2 /n term is absent). The plot 
shows clear differences in some regions of space, thus un- 
derlining the significance of the j 2 jn term. 

In our second example we send a fast, but still non- 
relativistic, proton against one of the carbons of ethene 
(Fig. [2J). The initial configuration is shown in panel (a). 
While approaching the carbon atom the proton accumu- 
lates some charge around it [Fig. Efb)]. It then scatters 
and leaves the system picking up part of the electronic 
charge. The electron-nuclei system is thus excited (in to- 
tal the electronic system absorbs around 20 eV). In panels 
(d) and (e) the leftmost carbon has already broken the 
two bonds with the hydrogens (that will later form an H2 
molecule). Clearly visible is also the polarization of the 
carbon-carbon double bond, and the starting of the for- 
mation of a lone pair above the leftmost carbon. We em- 
phasize once again that the formation of lone pairs can- 




FIG. 3: (Color) Snapshot of the time-dependent ELF with 
(upper panel) and without (lower panel) the j 2 /n term for 
the excitation of acetylene by a 17.15 eV laser pulse. The 
legend is the same as in Fig. Q] 
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not be visualized in movies of the time-dependent density. 
Only the TDELF allows one to observe this feature. At 
the end of the simulation [panel (f)] we can observe an 
H2 molecule (left), and two CH fragments (middle and 
right). The rightmost CH fragment is again breaking to 
yield a carbon and a hydrogen atom. Note again the 
lone pairs characteristic of CH, localized near the carbon 
atoms. 

These two examples illustrate the amount of informa- 
tion that can immediately be grasped just by looking at 
the time-dependent ELF, from 7T-7T* transitions, to the 
creation of lone pairs. One can infer the time scale and 
the temporal order of the various processes occurring, 
e.g., in a molecular collision: One can tell which bond 
breaks first and which breaks second, and how many at- 
toseconds or femtoseconds it takes to form new bonds and 
new lone pairs. We emphasize once more that with the 
advent of attosecond pulses this information will soon 



become experimentally available. We expect the time- 
dependent ELF to be a valuable tool in the analysis of 
many physical processes. One example is the creation 
and decay of collective excitations or the scattering of 
electrons from atoms and molecules. Another example is 
the process of vision: light promotes the electrons of reti- 
nal into the first excited state, which by its turn induces 
the isomerization of the photo-receptor and eventually 
leads to the firing of a neuron. The key feature of the 
time-dependent ELF, in all cases, is the time-resolved ob- 
servation of the formation, the modulation, or the break- 
ing of chemical bonds, thus providing a visual under- 
standing of the dynamics of excited electrons. 
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